Transition redshift from the V-reconstruction method 
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We show the effectiveness of using the procedure of the potential reconstruction (the V-method) 
in the estimation of the transition redshift corresponding to switch from deceleration to acceleration 
phase of the Universe. We investigate the FRW models with dark energy by applying the particle- 
like description of their dynamics. In this picture the evolution of the universe is represented by 
motion of a fictitious particle in the one-dimensional potential V{a) or V{z) where a is the scale 
■ factor and z is the redshift. The V-method solves the inverse problem where we reconstruct the 

' potential function from empirical data. We use the Riess et al.'s gold and gold-|-silver samples of SN 

, la data in this reconstruction. In the same framework we obtain both the estimates of the form of 

04 ■ the equation of state for dark energy and the transition redshift. We obtain that transition redshift 

5_( ' is O.SSlo Jy (gold subset) and 0.371q o7 (gold-|-silver sample) when the linear model of dark energy 

^ ll ('"^(2) ~ Wo -f wiz) and Qm,o ~ 0.3 are assumed. We compare the estimation of transition redshift 

. with the Riess et al. results. The V-method also alows us to find the value of Hubble function at 

the moment of transition. 
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I. INTRODUCTION 



. The recent supernovae observations indicate that our universe is currently accelerating QQ. However, from the 
' standard cosmological model we know that earlier the universe decelerated because of domination of matter, and 
. there should exist the moment when the switch from the deceleration phase to the acceleration phase takes place. 
The value of redshift for this transition zt has been recently found by Riess et al. They obtained zt = 0.46 ±0.13 
using a kinematic model independent on the content of the universe. Their method does not allow to obtain both the 
equation of state parameters and the transition redshift in the same framework. The dark energy equation of state 
O I parameters are estimated in models with different ansatz on w(z). However the constraint on the transition redshift 
^ ■ was obtained using other relation for dL{z). Let us note that there is no reason to treat the value of zt estimated in 
23 ' a certain model as a real value in another model. 

In this paper we propose the method which allows to calculate both equation of state parameters and the transition 
redshift in the same model. We estimate the value of zt through the potential function l/[a(z)] of the Hamiltonian 
system determining the evolution of the universe. In this approach the evolution of the universe is reduced to the 
5J] ' motion of a fictitious particle with a unit mass in a one-dimensional potential. The information about the influence of 
matter content on the dynamics of the universe is derived from the shape of the potential function y[a(z)]. Namely 
the scale factor a accelerates (decelerates) in the interval of z if I/[a(z)] is an increasing (decreasing) function of z. 
Therefore we should expect the potential function has a maximum for the universe with deceleration during its earlier 
epoch and acceleration at present. 

In our previous papers 0, Q we developed the potential function method to probe the dark energy. The use of 
the potential l/[a(z)] instead of wxiz) has an advantage because the former suffers less from the smearing effect 
caused by the double integral which relates wx{z) and ^^(z) 0. In these papers we showed that the dynamics of the 
Friedmann-Robertson- Walker (FRW) model of the universe, filled with non-interacting, non-relativistic matter with 
equation of state p„i ~ Wm{a)pm and dark energy with the equation of state px = wx{o,)px, can be reduced to the 
Hamiltonian system 

n^^ + Via), (1) 
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where an overdot denotes the differentiation with respect to the cosmic time t, and the potential 



V{a) 
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Poff(a) = Pm{a) + px{a), 



(2) 



where the units SttG = c = 1 are used and poff is the effective energy density of the mixture of non-interacting fluids 
wliich satisfy the conservation equation 



-i-{pi +Pi). 



(3) 



Then poff also satisfies equation (PJ. On the other hand the potential function V{a) can be reconstructed from the 
SN la data through the Hubble function H{z) which in the spatially flat universe is related to the luminosity distance 
dL{z) as 



H{z) 



d 

dz 



dL[z) 



(4) 



Because the trajectories of system always lie on the zero-energy level (as a consequence of the Hamiltonian 
constraint) we can rewrite potential to the form 



V{z) 



dz 



l + z 



(5) 



The potential V{z) contains all information necessary to determine the full dynamics of the universe on the phase 
plane (a, ct). Thus both the effective energy density Poff(a) and the equation of state coefficient Wcff[a(z)] = Pce/PcS 
can be calculated unambiguously from the potential function V{a) 



Poff(a) = - 



6Via) 



(6) 



where /y(a) = (dlnV) /{d\na) is the elasticity of the potential function V with respect to the scale factor a. 

As it is well known the behavior of the potential function in the neighborhood of a maximum can be approximated 
by a quadratic part of the expansion in Taylor's series. Therefore we have 



1 d'^V 



2 9a2 

where V{aT) — — dy/2 from the Hamiltonian constraint, and 



(a - arf 



Poff(a) = 



1 d^V 



2 da^ 



(7) 
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Pcs{z) = - 3(1 + zt)' 



+ -i{l + zf 



dz"^ 



(1 



ZT) 



{1 + zt)'' 



6(l + z)(l-h2T)^ 



dz^ 



(9) 



In the case of transition from deceleration to acceleration one can see from formula Q that there are three terms 
which can be interpreted as the cosmological constant (positive), the curvature contribution (negative), and the 
topological deffects with positive energy {p = —2/3/9). Note that as z is close to —1 (a ^ uq) than the last two terms 
are negligible and only the "cosmological constant" dominates the dynamics of the late epoch of the universe. 

The value of the potential function at the maximum is given by 



V{aT) = -lH^{l + ZT)-^. 



(10) 



The value of z^ and V{zt) inform us about the moment of transition and the Hubble function at this moment. The 
estimation of these parameters will be done in the next section. 
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II. ESTIMATION OF REDSHIFT TRANSITION 



In our previous papers the method of reconstruction of the potential function was applied to the analysis of global 
dynamics of the FRW models in the phase space 0, 0, ■ The potential function was reconstructed from the SN la 
data which allow us to determine all evolutional paths of the model for all admissible initial conditions. These results 
have the qualitative character, but it is also possible to obtain some quantitative attributes of the reconstructed 
dynamics of a model. 

It seems to be important to determine when the switch from the deceleration to the acceleration happens in different 
cosmological models. Because it is possible to differentiate between cosmological models through the quantitative 
features of their potential functions, we propose to use it as the test of cosmological models. The maximum of the 
potential function, interpreted as the time of switch between deceleration and acceleration phases, can be chosen as 
a one of the criteria (another may be the value of the potential function at zt). If someone obtains the value of zt 
independently from any model (e.g. one can numerically probe it directly from the observational data), it could be 
compared with values of zt found in different cosmological models. 

We calculate the probability distribution of the transition redshift for the spatially flat FRW model filled with 
nonrelativistic matter (p™ = Pm.oo-^^) a-nd dark energy for which the linear parametrization of the equation of state 
is assumed {w{z) = wa + wiz). In our analysis we use the recently released SN la data by Riess et al. Q. We use 
both the gold subset and gold+silver (full) sample. To estimate the cosmological model parameters wq, wi, and the 
transition redshift zt {^m,o = 0.3 and Hq = 65km s^^ Mpc^^ priors are used) we minimize the statistics given by 
the formula 



E 
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[fJ-t,jizj; ilm,0, Ho, Wo, Wi) - /ioj] 



(11) 



where fioj = Tricj ~ M is the extinction-corrected distance modulus for SN la at redshift Zj, cr^^j is the 
uncertainty in the distance moduli /io j including the dispersion in galaxy redhift due to peculiar velocities, 
fit,j{zj;n„i^o,Ho,wo,wi) = 5logdL{zj]n„ifi,Ho,wo,wi) + 25, 



dLiZj;Qrn.fi,Ho,Wo,Wi) = (1 + Zj) 

= (1 + z,)H^' 



dz' 



H{z';Q,n,o,Ho,wo,wi) 
dz' 



Vf7„,,o(l + z')^ + f^x,o(l + z')3(i+"'o-'"i) exp {3wiz) 



(12) 



The sumation in equation Ullf) is over all of the observed supernovae. The potential function for the model assumes 
the form 



V[a{z),n,n,o,Ho,wo,wi] = -^H'^{a{z),n,nfi, Ho,wo,wi)a'^{z) = -i 



d_ f dLjajz), ^},nfi, Ho,Wo,Wi) 
dz \ 1 + z 



a\z) 



1 



+ z) + nx.oil + z)3(i+-»— exp {iwiz) 



(13) 



TABLE I; Fitting results of the statistical analysis from distant type la supernovae data. 
(Sample of 186 SN la, Riess et al. (2004)) 



Sample 






wo 


Ml 




ZT 


V[zt) 


Method 


gold 


0.30 
0.30 


0.70 
0.70 


-1.35 

1 OO+0.19 


1.52 

1 40+0.81 


175.57 


0.40 

0.38l-.i? 


-1.60 

1 co+0.05 
1.00_Q 06 


best fit 
max(P) 


gold + silver 


0.30 
0.30 


0.70 
0.70 


-1.41 

^•■-'='-0.18 


1.73 

'U-O.SS 


229.36 


0.38 

Q 07 + 0.07 

'J-'J' -0.06 


-1.58 

1 rfi + 0.05 

1.00_o 04 


best fit 
max(P) 



The best fitted parameters of the model as well as the transition redshift zt and the maximum of the potential 
Vt = V{zt) are presented in Tabled 
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FIG. 1: The left column refers to the gold subset of Riess et al.'s SN la sample and the right column to the gold+silver sample. 
It is assumed that Qm,o = 0.30 and Ho = 65 km s~^ Mpc. Figures (a) and (b) present the levels of constant transition redshift 
on the plane (u;o,wi) as well as the confidence levels (68%, 95% and 99%) for pairs of {wq, wi) parameters obtained from the 
maximum likelihood method. The star denotes the best fit values for wo = —1.35, wi = 1.52 (gold subset) and for wo = —1.41, 
wi = 1.73 (gold+silver sample). The dot denotes the cosmological constant model. The shaded region is the la confidence 
level for zt from the likelihood method around the best fit value (the dashed line). Figures (c) and (d) present the probability 
distribution of the transition redshift zt- Figures (e) and (f) show the confidence levels for the reconstructed potential function 
V{z) as well as for the maximum of the potential function (zt, V{zt))- 
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FIG. 2: We draw the levels of constant zt on the plane (too,'"'i)- We distiguish zt = 0.46 (the dashed line) with its la 
confidence level. We also mark with the star the best fit of (i«o,wi) and the 68%, 95% and 99% confidence levels obtained by 
Riess et al. The key observation is that the transition redshift for their {wo,wi) parameters is lower than the corresponding 
value obtained from the Turner- Riess kinematic description |^. 



Assuming that the measurement uncertainties are Gaussian we can simply calculate the likelihood function 
^{^m,o, Hq, wq, Wi) from a chi-squared statistics in the following way 

rfn u \ ^ f x'^{^7n,o,Ho,wo,wi) 
C(ilmfi,Ho,wo,wi) = -^exp I 

N^JJJfc^p (- ^'^""'■°'f'"°'"^^ )d»^,orfgoc^^orf^i, (14) 

where N is the normalization coefficient given as an integral of the likelihood function over all probability variables. 
The confidence levels (68%, 95%, 99%) for pairs {wq.wi), obtained for the gold subset and full sample, are drawn in 
Fig. ni^a) and Fig. ^^b) respectively. We find the onc-dimensional probability distribution functions for wi and W2 
separately by integrating the likelihood function C over remaining probability variables and we obtain wo = — 1.33^q :[g, 

wi = lA9to.ll for the gold subset and wq = -1.39j:^:i^, ^"i = l-^Otass for t^e full sample. 

We find the probability distribution for the transition redshift zt by integrating the likelihood function C along 
the levels of constant zt (see Fig. ^c),(d)). As a result we obtain zt = O.SSIq'J!^ (gold subset) and zt = O.STlQQg 
(gold+silvcr sample). 

We also calculate the probability distribution for the value of the potential V{zt) by integrating the likelihood 
hmction C along the levels of constant V{zt) for the best-fitted value of zt- In this case we have V{zt) = -l-58to°ol 
in miits of IO^Hq (gold subset), V{zt) = — 1.58tg Q4 units of lO^H^ (gold+silvcr sample). 

Because of the Hamiltonian constraint we can also calculate the value of the Hubble parameter corresponding to 
transition redshift zt 



HT = il + ZT)V-^ViZT) (15) 

and the velocity of the scale factor during the transition epoch 

dr - V-2V{zt) (16) 

and obtain that Hx — 79.2 > Hq and dr — 56.6 < dg = Hq (gold subset). 

The quality of the fitting of the potential function as well as the position of its maximum is illustrated in Fig. ^ 
(e), (f). It is interesting that the value of V{zt) is calculated with lower error than zt- We also see that the errors 
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decrease in the neighbourhood of the region around the maximum of V{zt)- For the gold+silvcr sample the errors 
are lower than for the gold subset. 

We estimated the value of zt in the model with linear form of the equation of state coefficient w{z). The value of 
transition redshift was also obtained by Riess et al. in the different model. To compare their result in context of the 
presented framework we look for the transition redshift for their estimated model (model with linear ansatz on w{z)). 
We obtained lower value (marked as a star in Fig.[2Il than their corresponding value 0.46 (marked as the dashed line). 
In our opinion the reason of this difference is that their estimation of zt was done in different model without any 
priors on model parameters. 

III. CONCLUSIONS 

We presented the V-method of reconstruction of the potential function. We applied it to the estimation of the state 
parameters of the cosmological model with dark energy. We concentrated on the estimation of the position of the 
maximum of the potential function which can be interpretated as the moment of switch from the deceleration to the 
acceleration phase during the Universe evolution. 

The shape of the potential function is reconstructed from the distant SN la data. For our analysis we used the 
samples of SN la prepared by Riess et al. 

The Hamiltonian formulation of the dynamics enables us to determine whole dynamics of the cosmological model 
from the potential function only. It is interesting that this function determines not only the qualitative structure of 
the phase space 0, but also gives us some quantitative information about the model parameters. 

We showed that the V-metliod is very effective in analysis of cosmological models with dark energy where we can 
find the moment of transition as the quantitave parameter. The other obtained parameter of state of the system is the 
velocity of the scale factor (or the Hubble parameter) at the moment of transition. We can treat these two estimated 
parameters as initial conditions and determine uniquely the evolution path of the model. 

The main results of using this method to the cosmological model with dark energy are 

1. The method allows to study quantitative and qualitative aspects of dynamics of models. 

2. We showed that the state parameters at the moment of transition can be obtained in any models for which the 
reconstruction can be done. 

3. We estimated the value of the transition redshift and the value of Hubble function at the transition moment 
moment. 

Riess et al. estimated the model parameters assuming the different forms of the equation of state coefficient w{z), 
but to find the transition redshift zt they used a certain kinematic model with a linear expansion for q{z). It seems 
to be no a priori reasons to extrapolate the value of zt between two different models. The question is to which w{z) 
model corresponds their value of zt- 

Our methods overcomes this obstacle. To obtain the model parameters and the transition redshift we use the same 
dL{z) relation. 
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